Non-Langevin behaviour of the uncompensated magnetization in nanoparticles of 

artificial ferritin. 
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The magnetic behaviour of nanoparticles of antiferromagnetic artificial ferritin has been inves- 
tigated by 57 Fe Mossbauer absorption spectroscopy and magnetization measurements, in the tem- 
perature range 2.5K-250K and with magnetic fields up to 7T. Samples containing nanoparticles 
with an average number of 57 Fe atoms ranging from 400 to 2500 were studied. By analysing the 
magnetic susceptibility and zero field Mossbauer data, the anisotropy energy per unit volume is 
found, in agreement with some of the earlier studies, to have a value typical for ferric oxides, i.e. a 
few 10 s ergs/cm 3 . By comparing the results of the two experimental methods at higher fields, we 
show that, contrary to what is currently assumed, the uncompensated magnetisation of the ferritin 
cores in the superparamagnetic regime does not follow a Langevin law. For magnetic fields below 
the spin-flop field, we propose an approximate law for the field and temperature variation of the 
uncompensated magnetization, which was early evoked by Neel but has so far never been applied 
to real antiferromagnetic systems. More generally, this approach should apply to randomly oriented 
antiferromagnetic nanoparticles systems with weak uncompensated moments. 
PACS: 75.50 Tt; 75.50 Ee; 76.80 +y 
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I. INTRODUCTION 

Natural ferritin is the iron-storage protein of animals, 
plants and bacteria. It is composed of a ferrihydrite-like 
core with formula (FeOOH^FeOH^PC^) about 7nm in 
diameter and containing up to about 4500 Fe atoms |Q, 
surrounded by a 12 nm diameter multisubunit protein 
shell. By a suitable chemical synthesis process, it is pos- 
sible to reconstitute the ferrihydrite ferritin core inside 
the empty protein shell (in this case obtained from horse 
spleen ferritin by reductive dissolution of the native core) 
. This way, one can monitor the amount of iron avail- 
able for the build up of the core, and a better control 
of the average core size is obtained. A number of in- 
vestigations have been performed, both on natural and 
artificial ferritin, using either 57 Fe Mossbauer absorp- 
tion spectroscopy or magnetization measurements 
Jl-El). They have shown that the magnetic structure of 
the Fe 3+ ions in the core is probably antiferromagnetic, 
and that the particles possess a (small) uncompensated 
magnetic moment. However, the interpretation of the 
data obtained by the two techniques is somewhat contra- 
dictory. Indeed, on one side, the Mossbauer studies with 
applied magnetic field reveal that, for randomly ori- 
ented ferritin samples, the Fe 3+ magnetic moments in the 
cores remain close to their easy axes, even in a strong 
field of 14 T. On the other side, the magnetic measure- 
ments have been interpreted so far by assuming that the 
uncompensated magnetization follows a Langevin law, 
which implicitly assumes the Fe 3+ moments are free to 
rotate. Besides its fundamental interest, the understand- 
ing of the magnetic properties of ferritin is of importance 
for NMR relaxometry [|o| and biophysical purposes as it 



is a good contrasting agent for NMR imaging. The fer- 
ritin core has also recently arisen great interest because it 
has been thought to be a good candidate for the observa- 
tion of quantum tunneling of the magnetization p2[ , and 
various experimental studies, at low temperature, have 
claimed to have observed this tunneling 0-na|. 



The aim of the present work is to reexamine the prob- 
lem of the uncompensated magnetization in ferritin, by 
comparing 57 Fe Mossbauer data, both in zero field and 
in applied fields up to 7 T, with magnetization data ob- 
tained for the same samples. We show that, by taking 
proper account of the crystalline anisotropy in an anti- 
ferromagnetic structure, one obtains a coherent interpre- 
tation of the results of the two techniques. We present a 
new interpretation of the superparamagnetic behaviour 
of an ensemble of randomly oriented antiferromagnetic 
particles, in terms of a non-Langevin law for the uncom- 
pensated magnetisation, for applied magnetic fields lower 
than the spin-flop field. This approach, which was evoked 
by Neel in an early work ]lq| , has however never been 
considered in the modern and quantitative magnetic in- 
vestigations of ferritin. Another motivation of our work 
is to examine the influence of the mean iron content of 
the (artificial) ferritin core (ranging from 400 to 2500 Fe 
atoms) on the physical properties, such as the anisotropy 
energy, or the Neel temperature. 

The paper is organized as follows: section [n] contains 
the description of the artificial ferritin samples, section 
III briefly recalls the superparamagnetic behaviour of an- 
tiferromagnetic particles, both from the point of view 
of magnetometry and 57 Fe Mossbauer spectroscopy, sec- 
tions IV and [y| contain respectively the FC-ZFC mag- 
netic susceptibility and zero field Mossbauer measure- 
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merits, and the high field Mossbauer data; in section 



VI, the influence of anisotropy in uncompensated anti- 



ferromagnets is examined, and section VII reports the 



magnetisation experiments and the new interpretation 
we propose. 



II. SAMPLE CHARACTERIZATION AND 
EXPERIMENTAL TECHNIQUES 

In order to investigate the effect of decreasing both 
the particle size and the mean number of Fe atoms per 
ferritin core, the experiments were performed with a set 
of artificial ferritin samples with different Fe loadings. 
Five samples were studied, having an average iron load- 
ing ranging from 400 to 2500 Fe ions per ferritin core. 
In order to obtain a good signal to noise ratio in the 
Mossbauer spectra, the Mohr salt used as the starting 
material was prepared with iron 95% enriched in 57 Fe. 
The protein concentration was determined by the Lowry 
method and the mean iron content by atomic adsorption 
analysis. In the following, the samples will be labelled 
according to the mean number of Fe atoms per core de- 
termined by this latter method. The samples consist of a 
solution of artificial ferritin with a concentration of about 
lOmg/ml, that is ten times smaller than the concentra- 
tion of commercial ferritin (about lOOmg/ml). Then, 
the distance between particles in the solution is about 
50 nm, which corresponds to a dipolar field of magnitude 
10 mG (at saturation of the uncompensated moments). 
Our samples can therefore be considered as an ensemble 
of non-interacting particles. Pictures taken by transmis- 
sion electron microscopy (TEM) show the cores to be 
discrete and roughly spherical in shape. Size histograms 
have been established by measuring the diameter of 500 
particles taken from different parts of the grid. For each 
sample, the diameter distribution is rather narrow and 
can be fitted to a lognormal shape, with a mean diameter 
do ranging from 4 nm for the particles with 400 Fe atoms 
per core to 5.7 nm for the particles with 2571 Fe atoms 
per core, and a standard diameter deviation a about 0.15. 
Electron diffraction patterns, performed on the particles 
with 2571 Fe atoms per core, give d-spacings correspond- 
ing to well-ordered ferrihydrite. 

The magnetisation measurements were made with a 
commercial SQUID magnetometer in magnetic fields up 
to 5.5 T in the temperature range 2.5K-250K. For all 
temperatures and fields, we have measured both the sig- 
nal of the solution containing ferritin and the signal of 
the solution containing apofcrritin (the empty protein 
shells) with the same concentration. After subtraction 
of the second signal from the first, we thus obtain the 
signal due only to the ferritin cores. For the low field 
(80 G) magnetisation measurements, both Field Cooled 
(FC) and Zero Field Cooled (ZFC) branches were mea- 
sured. In the Zero Field Cooled (ZFC) procedure, the 
sample was cooled in zero field from room temperature 
down to 2.5 K, and for the Field Cooled (FC) branch, the 



sample was cooled with a field of 80 G from room temper- 
ature to 2.5 K; in both cases, measurements proceeded on 
heating. 

The 57 Fe Mossbauer absorption spectroscopy experi- 
ments were performed both in zero magnetic field and 
in magnetic fields up to 7T applied perpendicular to 
the 7-ray propagation direction, between 4.2 K and 90 K. 
The Mossbauer spectra were obtained using a 57 Co:Rh 
source, mounted on an electromagnetic drive with a tri- 
angular velocity signal. The ferritin solution was placed 
in a copper holder covered with pure thin Aluminium 
sheets, ensuring a good thermalization. For the in-field 
measurements, the ferritin solution is first frozen in zero 
field and therefore there is no preferential orientation of 
the magnetization of the particles. 



III. SUPERPARAMAGNETISM OF 
ANTIFERROMAGNETIC PARTICLES 



The Fe 3+ ion has a saturated magnetic moment too 
=5 fiB- The magnetic structure of the ferritin cores is ex- 
pected to be antiferromagnetic, with a Neel temperature 
Tjy of the order of a few hundred Kelvins. In zero external 
field, the preferential orientation of the two antiferromag- 
netic sublattices in the core is determined by the crys- 
talline anisotropy ("antiferromagnetic axis"). For anti- 
ferromagnetic nanometric particles, one expects a small 
uncompensated magnetization which may arise from the 
core (presence of defects) and/or from the unpaired sur- 
face moments of the particle ||l7| . Due to the strong 
exchange interaction, the uncompensated moments are 
aligned with the antiferromagnetic sublattice moments, 
except probably at the surface due to broken exchange 
bonds (lq . At a given temperature T below Tn, the 
moments of the two sublattices fluctuate by crossing the 
anisotropy energy barrier: this is the superparamagnetic 
relaxation. In the case of an axial anisotropy, the relax- 
ation time for the reversal of the direction of the magne- 
tization of a particle with volume V is described by an 
Arrhenius type equation first proposed by Neel |19|: 



KV 
T ° eXP< -fc^T' ) ' 



(1) 



where K is the magnetic anisotropy energy per unit vol- 
ume, To a microscopic relaxation time usually considered 
to be a constant with magnitude 10 -10 -10~ n s but vary- 
ing in fact with V and T [EJ. Due to the particle size 
distribution present in any real sample, there is a dis- 
tribution of anisotropy barriers KV which, according to 
Eqnjl], results in a broad distribution of relaxation times. 
For a given measurement technique, with a characteristic 
time T m , the blocking volume: 



k B T 
K 



ln(r m /r ), 



(2) 
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defines two populations of particles: those with V > V& 
and those with V < V&, whose magnetization fluctuation 
time is respectively larger and smaller than r m . Using 
Eqn.(Q), one can also define a mean blocking tempera- 
ture: 



3(1 + cos 2 4>) : 4 sin 2 </>:(!+ cos 2 



(6) 



K(V) 



k B ln(T m /r ) : 



(3) 



such that, for T <C Tb all the particles are in the "frozen 
regime" and for T ^> Tf, all the particles are in the "su- 
perparamagnetic regime" . 

For magnetization experiments, the characteristic time 
can be considered to be the measurement time t x , which 
is about 100 s. In low fields, and for the superparamag- 
netic regime, Neel has shown that the magnetization is 
the sum of two contributions [Elf : 



M(H,T) = X A F H + M nc (H,T), 



(4) 



where xaf is the susceptibility arising from the canting 
of the two antiferromagnetic sublattices, and M nc is the 
magnetization due to the uncompensated moments. At 
T=0 and in the case the magnetic field is perpendicular to 
the antiferromagnetic axis, xaf = X± which, in the stan- 
dard molecular field theory, is given by: x± — Mo/ He,, 
Mq being the magnetization of a sublattice and H-e the 
exchange field. In the low field limit, the uncompen- 
sated magnetization follows a Curie law, irrespective of 
the strengh of the anisotropy pi] . At higher fields, the 
choice of an expression for M nc (H,T) is not a trivial 
problem, and an important issue of the present work is 
the determination of the best approximation for the field 
and temperature dependence of the uncompensated mag- 
netization in antiferromagnetic particles. 

The characteristic Mossbauer time is the Larmor pe- 
riod tl associated with the magnetic hyperfine interac- 
tion. For 57 Fe, it is of the order of 5xl0~ 9 s. In the case 
of magnetic ordering of the Fe 3+ ions, there is a magnetic 
hyperfine field Hhf at the nucleus, proportional to the 
Fe 3+ moment and directed opposite to it. If the fluctua- 
tion time r of the magnetization (and hence of the hyper- 
fine field) is longer than tl, the Mossbauer spectrum is 
a six- line hyperfine field pattern; if it is smaller than tl, 
the Mossbauer spectrum is a two-line quadrupolar pat- 
tern. In real systems, due to the broad distribution of r 
values, a zero field pattern consists of a superposition of 
magnetic and quadrupolar subspectra, in a sizeable tem- 
perature range. When a magnetic field is applied, each 
nucleus experiences an effective field H e ff , which is the 
vectorial sum of the applied field H and of the hyperfine 
field H hf : 



H e ff — H + H 



hf • 



(5) 



If <fi is the angle between the 7-ray propagation direction 
and the effective field direction, the intensity ratios of the 
outer to middle to inner pairs of lines of the Mossbauer 
sextet are given by: 



If the effective field is aligned perpendicular to the 7-ray 
direction, the intensity ratios are 3:4:1, and for a random 
orientation of the effective field, the intensity ratios are 
3:2:1. 



IV. LOW FIELD SUSCEPTIBILITY AND ZERO 
FIELD MOSSBAUER EXPERIMENTS 



The two FC and ZFC branches of the magnetic sus- 
ceptibility for the particles with 410 Fe atoms per core, 
measured in a field ff=80G, are represented in Figj^. 
The signal is almost entirely due to the uncompensated 
moments, the antiferromagnetic s usce ptibility being very 



weak as will be shown in section VII 
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FIG. 1. Thermal variation of the FC and ZFC magnetic 
susceptibility in the artificial ferritin sample with a mean Fe 
loading of 410 atomes per core, measured in a field of 80 G. 
The solid lines are fits using to Eqn.fy and a Gaussian dis- 
tribution of u, nc values. 

The shape of the curves is typical of an ensemble of 
non-interacting relaxing moments with a distribution of 
anisotropy barriers p2fl . The peak temperature of the 
ZFC curve for the particles with 410 Fe atoms per core 
is Tpeak ~7.5K, and the irreversibility temperature at 
which the FC and ZFC branches join is Ti rr ~ 15 K. To 
interpret these curves, we follow here the model of Rcf. 
|p2| . Assuming that, for a given particle volume, the un- 
compensated moment iL nc is a known function of V, the 
thermal variation of the FC and ZFC magnetic moment 
writes: 



m nc (T) 



H 



V b (T) 



3k B T j Vr 
H ' v < 



n. nc (V) 2 f{V)dV+ 



"iK 



V b (T) 



(7) 
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where /(V) is the log-normal volume distribution func- 
tion, Vb(T) the blocking volume associated with t x , and 
with To = 1CP 11 s, and 77=0 for the ZFC branch and In ^ 
for the FC branch. However, in antiferromagnetic par- 
ticles, the relationship between volume and uncompen- 
sated moment is not well known. The most probable is 
that, for a given volume, there is a distribution of uncom- 
pensated moment values. According to an argument due 
to Neel Q , the mean value of this distribution is propor- 
tional to N p , where N is the number of Fe atoms in the 
particle and p=l/2, 1/3 or 2/3 in the case respectively of 
a random distribution of the uncompensated moments in 
the volume, at the surface, or in the presence of regular 
"active" (ferromagnetic) planes at the surface. It appears 
that the fit of the experimental FC and ZFC curves us- 
ing Eqn.g is not very sensitive to the particular shape of 
the distribution of uncompensated moments, nor to the 
choice of a particular p value. The solid lines in Fig|l], 
which reproduce quite well the experimental data, were 
obtained by introducing a Gaussian distribution of ii nc 
values into expression (R). 




FIG. 2. Thermal variation of the magnetization of artifi- 
cial ferritin samples with different Fe loadings (the figure near 
each curve represents the mean number ofFe atoms per parti- 
cle) upon warming in a field of 50 G after cooling in zero held 
(ZFC curves). The fit using Eqn.^) (solid lines) is shown 
only for the sample with a loading of 2571 Fe atoms per core. 
The magnetization is expressed in emu per mg of protein. 



The mean value of this distribution is taken as: (fJ, nc ) — 
moiiy/N (assuming the Fe atomic density n = N/V is 
constant for a given particle set), and its standard devi- 
ation as a nc =0A (/x nc ). The distribution is restricted to 
positive fi nc values. The fitted value of the proportion- 
ality parameter u is in the range 0.8-1 for the different 
particle sets, and the standard deviation of the diame- 
ter distribution was kept at the value cr=0.15 obtained 
from the size histograms. This choice of a N 5 law for the 
mean fi nc value corresponds to a volumic random distri- 
bution of the Fe atoms bearing the uncompensated mo- 



ments but, as emphasized before, the data do not allow 
us to exclude the presence of uncompensated moments at 
the surface. Nevertheless, these fits allow the anisotropy 
constant K to be determined: we find that K lies in 
the range 3-6 10 5 ergs/cm 3 for the different particle sets, 
with a tendency to increase as the particle size decreases. 

In Fig]| IS re P resen ted, for each sample, the behaviour 
of the ZFC magnetization as a function of temperature 
(the fit is shown as a solid line only for the particles 
with 2571 Fe atoms per core). The peak temperature 
increases with the particle size, and it varies approxi- 
mately linearly with the mean barrier energy K(V), as 
can be seen in Fig||. As, in antiferromagnetic particles, 
the shape anisotropy is very small , this linear correla- 
tion is expected if the volumic anisotropy dominates over 
the surface anisotropy. Indeed, in this case, T pea k can be 
shown to be proportionnal to the blocking temperature 
Tb Plpl: 



I^peak — 9( a )Tb, 



(8) 



where g(a) is an increasing function of the standard di- 
ameter deviation a with g(a — 0) = 1. As a is relatively 
constant for the different particle sets, the peak tempera- 
ture must be proportional to the mean anisotropy barrier 
owing to Eqn.||. 
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FIG. 3. Dependence of Tp Ca it of the ZFC curves on the 
mean anisotropy energy in artificial ferritin samples. 

For all samples, the FC curve joins the ZFC curve 
above an irreversibility temperature which is about twice 
the peak temperature; therefore, for all particle sizes, the 
superparamagnetic regime associated with the magneti- 
zation measurement time occurs above 30 K. 

Representative 57 Fe Mossbauer absorption spectra in 
zero field for the particles with 410 Fe atoms per core 
are shown in FigJ^. At 4.2 K, a static magnetic hyperfine 
pattern is observed with a hyperfine field Hhf — 490 kOe, 
identical for all particle sets. As temperature increases, 
the quadrupolar doublet progressively replaces the six- 
line magnetic pattern as the superparamagnetic fluctu- 
ation time of more and more particles becomes smaller 
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than tl, i.e. as the blocking volume given by Eqn.(g) as- 
sociated with the characteristic time tl becomes smaller. 
The spectra were fitted to a superposition of a quadrupo- 
lar hyperfine spectrum and of magnetic hyperfine spectra 
with a distribution of hyperfine fields (histogram) . 
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FIG. 4. 5 Fe Mossbauer absorption spectra at selected tem- 
peratures in zero Geld in the artificial ferritin sample with a 
mean Fe loading of 410 atoms per core. 

The thermal variation of the relative intensity f p (T) of 
the quadrupolar doublet derived from these fits is shown 
in Fig.|^. This "paramagnetic fraction" can be calculated 
in the frame of Neel's model of thermally activated fluc- 
tuations: 



1 fV b (T) 

f P (T) = W} J Vf(V)dV, 



(9) 



where V b (T) = In ^ and t = 10" n s. The para- 
magnetic fraction is independent of the values (or distri- 
bution) of the uncompensated moments as long as one 
neglects the possible (but unknown in antiferromagnets) 
dependence of To on /i„ c . For all particle sets, the ther- 
mal variation of f p is well reproduced using Eqn.^ (solid 
line in Fig.^), with an anisotropy constant close to that 
derived from the fit of the FC-ZFC susceptibility curves. 
The blocking temperature, defined as the temperature for 



which half of the particles are in the superparamagnetic 
regime, increases with the particle mean size. For all iron 
loadings, all the particles are in the superparamagnetic 
regime above 80-100 K. 
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FIG. 5. Thermal variation of the relative intensity of the 
quadrupolar doublet in the artificial ferritin sample with a 
mean Fe loading of 410 atoms per core. The solid line is a fit 
using Eqn.(Q). The blocking temperature is 22 K. 



V. HIGH-FIELD MOSSBAUER EXPERIMENTS 

According to the zero field Mossbauer data, for all par- 
ticle sizes, the fluctuations of the magnetization at 4.2 K 
are slow with respect to the hyperfine Larmor frequency 
and they are fast at 90 K. In-field Mossbauer spectra at 
these temperatures have been recorded in natural fer- 
ritin and haemosiderin in Ref. ||. We repeated these 
measurements in our artificial ferritin samples to try and 
detect the influence of the particle size on the orientation 
of the magnetic moments in the frozen and paramagnetic 
regimes when a magnetic field up to 7 T is applied. We 
actually found that there is no drastic qualitative change 
of the in-field spectra for the different particle sets, and 
our data are similar to those recorded in natural ferritin. 
We therefore only present the spectra for the particles 
with 410 Fe atoms per core, and our interpretation fol- 
lows the same lines as that made in Ref. 



A. Spectra in the frozen regime (T=4.2 K) 

The spectra at 4.2 K for fields of 3 T and 7 T in these 
particles are shown in Fig.[| They are resolved six line 
hyperfine spectra, analogous to the zero-field spectrum 
(see top of FigJ^). The intensity ratios of the pairs of 
lines remain close to 3:2:1, even at 7T, which means that 
the field does not drastically modify the orientations of 
the individual Fe 3+ moments with respect to those in 
the randomly oriented zero-field configuration. The spec- 
tra were then fitted assuming random orientations of the 
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static hyperfine field Hhf, with fixed magnitude. For a 
given angle 9 b between Hhf and the applied field, the 
magnitude of the effective field H e ff = H + Hhf experi- 
enced by the nucleus is given by: 



H eff = {(H hf cos 9 b + H) 



H 2 hf sm 2 9 b }i. 



(10) 



The orientational disorder of Hhf yields then a small dis- 
tribution of the magnitude of H e f f . It entirely accounts 
for the broadening of the lines with respect to the zero- 
field linewidths, which is the most visible at 7 T. 
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FIG. 6. Fe Mossbauer absorption spectra in the frozen 
regime (4.2 K) in the artificial ferritin sample with a mean Fe 
loading of 410 atoms per core, with magnetic fields of 3T and 
7 T applied perpendicular to the 7 -rays propagation direction. 
The solid lines are a fit to a model with random orientation 
of the hyperfine fields. 

As seen in Fig.|^, the model of a random orientation of 
the moments satisfactorily reproduces the shape of the 
Mossbauer spectra, except for a small misfit of the shape 
of the intermediate pair of lines. This implies that, up 
to a field of 7T, the magnetic moments of the Fe 3+ ions 
remain close to the antiferromagnetic axis in our ferritin 
samples at 4.2 K. The same qualitative conclusion can 
be drawn from the 14 T spectrum in horse spleen ferritin 
ofRef. §. In the following, we will label this moment 
configuration as the "random magnetic orientation" con- 
figuration. 



B. Spectra in the superparamagnetic regime 
(T=90K) 

In the superparamagnetic regime and in zero applied 
field, the fast fluctuations of the hyperfine field smear 
out the magnetic hyperfine structure. At 90 K, the zero 
field spectrum is a quadrupolar doublet identical with 
the main component of the spectrum at 35 K shown at 



the bottom of Fig.^. If a field is applied, the shape of the 
spectra, shown in FigJ^, changes, and at 5 T a six line 
hyperfine field pattern is distinguishable, but with large 
broadenings. In order to interpret the spectra, we will 
assume that a dynamic "random magnetic orientation" 
model holds at 90 K, i.e. that the moments fluctuate 
along or close to the antiferromagnetic axis. Our analy- 
sis follows here that of Ref. ||. For a given orientation 
9 b of the antiferromagnetic axis and a given value of the 
uncompensated moment, the energies of the two particle 
states with "up" and "down" orientations. of the uncom- 
pensated moment are separated by a Zeeman splitting: 



AE Z = 2fi nc (V)Hcos9 b . 



fill 



At a given Fe site, the hyperfine fields associated with 
each energy level are opposite, resulting in a thermally 
averaged hyperfine field: 



Hhf (T, V) — H^(T) tanh[ 



u, nc (V)H cos 9 b 



k B T 



(12) 



The values 



where H®JT) is the local hyperfine field, 
of Hh f {T, V) for a given particle set are therefore dis- 
tributed due to both the distribution of fJ. nc {V) values 
and to the random distribution in 9 b values. 




-5 5 
v (mm/s) 



FIG. 7. 57 Fe Mossbauer absorption spectra in the super- 
paramagnetic regime (T—90K), in the artificial ferritin sam- 
ple with a mean Fe loading of 410 atoms per core, with mag- 
netic fields of 2T and 5 T applied perpendicular to the 7-rays 
propagation direction. The solid lines correspond to fits us- 
ing the dynamic "random magnetic orientation" model (see 
text). 

The fits of the in-field spectra with this dynamic "ran- 
dom magnetic orientation" model, like those of the FC- 
ZFC curves, are rather insensitive to both the assump- 
tion made about the distribution of uncompensated mo- 
ment values and the value of the exponent p. For the 
fits shown as solid lines in Fig.^, we used the same dis- 
tribution of uncompensated moments as for the fits of 
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the FC-ZFC curves, i.e. a truncated gaussian shape with 
a mean value: (/i„ c ) = itiqu^/N. We find that the pa- 
rameter u is 0.4-0.5 for the lowest fields, and that it de- 
creases as the field increases. The good reproduction of 
the spectral shape, especially if one considers that u is 
the only adjustable parameter, indicates that it is real- 
istic to consider that the individual moments, and hence 
the sublattice magnetizations, fluctuate along or close 
to the antiferromagnetic axis at 90 K. The fact that the 
mean uncompensated moment found here is about twice 
smaller than that derived from the FC-ZFC curves prob- 
ably originates from the roughness of the assumption of a 
gaussian distribution of ii nc values. The decrease of u as 
the field increases, also observed in Ref. Q, is probably 
due to the fact that the "random magnetic orientation" 
model starts to break down at the higher fields, as will 
be discussed in the next section. The values obtained 
for u at the lowest fields (~0.5) are similar to those de- 
rived fr om t he magnetization measurements described in 
section 
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VI. INTERPLAY OF THE EXCHANGE, ZEEMAN 
AND ANISOTROPY ENERGIES IN 
UNCOMPENSATED ANTIFERROMAGNETS 

Before proceeding with the experimental results of the 
magnetization measurements and their interpretation, we 
will examine the effect of a magnetic field on an uncom- 
pensated antiferromagnetic particle, and determine the 
field range for which the "random magnetic orientation" 
configuration is expected to be a correct approximation. 
For a fully compensated antiferromagnet, the relevant 
threshold field is the so-called spin-flop field which writes: 
H s f ~ \/2HaHb, where He is the exchange field and 
Ha = the anisotropy field, Mq being the magnetiza- 
tion of one sublattice |Q . When the applied field H is 
aligned along the antiferromagnetic axis, the sublattice 
magnetic moments rotate to a direction perpendicular to 
it for H = H s f. For a particular orientation 6^ of the 
antiferromagnetic axis with respect to H, the T=0 equi- 
librium direction 9 of the sublattice moments, in the limit 
H/H E < 1 and H A / H E <C 1, is given by @: 



cos u = — 



cos 2 9 h 



^/(l + x 2 ) 2 -Ax 2 cos 2 9 b ' 



(13) 



where x = -§—■ For all initial orientations, a crossover to 
a "quasi spin-flop" configuration, where the equilibrium 
orientations are close to 90°, occurs for H > H s [ and, for 
lower fields, the equilibrium position 9 remains close to 
the initial orientation 9 b . 

In order to investigate the effect of the uncompensated 
moment, we will follow a model developed by Mjrup p7[ . 
For a particle with volume V and "degree of uncompen- 
sation" a = ^fy , the energy of an antiferromagnetic 
particle is given by: 



E(9, 6 b ) = M He[-(1 + a) + ^[3 2 sin 2 9+ 

af3cos9 + n(l+a)cos 2 {9-9 b )l (14) 



where j3 = H j He and k = Ha/ He- This expression is 
valid for (3 <C 1, k«1 and when the "canting angle" 2e 
of the two sublattices has its equilibrium (small) value: 
2e = j3 sin 9. Mjrup has shown that the spin-flop field in 
the presence of an uncompensated moment is enhanced 
with respect to H s f pTj: 



H% = —ocHe 



-a 2 H 2 
4 E 



(15) 



In our ferritin samples, using for instance the values 
for the sample with 410 Fe atoms per core: K ~ 
5 x 10 5 ergs/cm 3 and the atomic density: N ~ 1.2 x 
10 22 Fe 3+ /cm 3 , the anisotropy field is: H A ~0.18T. The 
exchange field is more difficult to estimate because vari- 
ous values for Tn have been given in the literature. We 
will take here our estimation of the Neel temperature, 
which we s how is close to 500 K for all particle sizes 
(see section |VIl| ) . Using the simple molecular field re- 
sult: fc^TN = ^g(S + l)i?EMSi one obtains an exchange 
field He — 320 T. As to the values of the degree of un- 
compensation a, they are probably distributed within a 
given particle set; using as a reasonable estimate of the 
mean uncompensated moment: ([i nc )= 0.5 moV^V, we 
find that the mean value of a ranges from 2% for the 
biggest particles to 5% for the smallest ones. 
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FIG. 8. Variation of the equilibrium orientation 9 of the 
sublattice magnetizations in an uncompensated antiferromag- 
net as a function of the orientation 9), of the antiferromagnetic 
axis with respect to the applied held, for He=320 T. The ra- 
tio of the uncompensated moment to the sublattice magnetic 
moment is: a=0.03. The curves (solid lines) are calculated 
for different values of the applied held H ranging from 3T to 
30 T. The dashed line is the curve 9 = 9b- 
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In Fig.|8j, we show the curves giving the T = equilib- 
rium orientations 9 as a function of the initial orienta- 
tion 9b, for the representative value a = 0.03, and for 
different values of the applied field H. For this a value, 
the spin-flop field is: #^=16.6 T. For fields below 7T, 
Fig.g shows that the equilibrium orientation 9 does not 
strongly depart from 9b for all initial orientations. Thus 
the T = "random magnetic orientation" model is seen 
to be a good approximation for fields lower than the max- 
imum field of our experiments (7T), thus justifying the 
use of this approximation to account for the T=4.2K in- 
field Mossbauer spectra. This model starts to break down 
for a field around 7T and, above the spin-flop field, the 
equilibrium 9 value is seen to be quasi independent of 9b 
and growing towards 90° as the field is further increased. 

For a fully compensated antiferromagnetic particle, 
the energy profiles present two equally deep potential 
wells, separated by 7T, corresponding to the invariance 
of the system by inversion of the two sublattice mo- 
ments. When the particle possesses an uncompensated 
moment, this symmetry is broken and the energy profiles 
become asymmetrical, as shown in Figj^ which represents 
E(9,9b = 50°) from Eqn.[l4|for different applied fields up 
to 15 T, and for a=0.03, i? E =320T and ff A =0.18T. The 
two potential wells are still clearly present for fields be- 
low 5 T, the energy difference between the lowest points 
in each well being essentially the Zeeman splitting given 
by Eqn.|Tl|. At higher fields (7 - 10 T), the less shallow 
well is somehow smeared out and it appears anew when 
the field is further increased. 
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FIG. 9. Energy profile for an uncompensated antiferromag- 
net, for 8b =50° , a degree of uncompensation ct=0.03 and an 
exchange Held He= 320 T, for different values of the applied 
held. The other parameters correspond to the particles with 
410 Fe atoms per core. 



ated with the canting of the two antiferromagnetic sub- 
lattices, \±H 2 - When a is not too small, this yields the 
approximate boundaries: Si < fj < aH^, i.e. this re- 
gion lies between 6 and 9.6 T for a =0.03. At low fields 
(H < -^fO, the positions of the two wells remain close to 
9b for the shallowest (see Fig]s|) and it — 9b for the other. 
Therefore, in this field region and in the superparam- 
agnetic regime, the dynamic "random magnetic orienta- 
tion" model, where the uncompensated moment is taken 
to fluctuate along the antiferromagnetic axis, is a fairly 
good approximation. At higher fields, the positions and 
shape of the potentiel wells start to be strongly modified, 
which probably accounts for the difficulty of coherently 
reproducing the Mossbauer spectra at 90 K using this 
model. 



VII. ISOTHERMAL MAGNETIZATION 
EXPERIMENTS 

Isothermal magnetization measurements have been 
performed in the superparamagnetic regime, between 
35 K and 250 K, for fields up to 5.5 T. Representative 
curves are shown in Figs.^ and for the particles with 
982 Fe atoms per core. For a given temperature, the 
magnetization increases with the field, without saturat- 
ing, and, for a given field, it decreases with increasing 
temperature, as observed in Refs. ||^— |Tl|] ■ 
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H (Tesla) 

FIG. 10. Low field part of the isothermal magnetization at 
35 K and 80 K in the artificial ferritin sample with a mean Fe 
loading of 982 atoms per core. The solid lines are fits with 
the "random magnetic orientation" model (see text), and the 
dashed lines represent the Bts with the Langevin model. The 
inset shows the G(x) function (see text) and the Langevin 
function L(x). 



The field range where only one potential well is present 
is such that the Zeeman energy associated with the un- 
compensated moment \x nc li is larger than the anisotropy 
energy KV , but smaller than the Zeeman energy associ- 



For these antiferromagnetic particles, the magnetization 
can be thought to arise from two contributions, as ex- 
plained in section III: a linear term xafH, where %af 
is the powder antiferromagnetic susceptibility and ac- 
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counts for the weak canting of the two sublattices, and 
another term due to the uncompensated moments. In the 
previous magnetization studies of ferritin quoted above, 
this latter contribution has been interpreted using the 
Langevin model. However, for such low fields, the dis- 
cussion presented in section VI shows that the uncom- 
pensated moments must be considered as fluctuating es- 
sentially along the antiferromagnetic axis, which is con- 
firmed by the analysis of the in-field Mossbauer spectra 
at 90 K. This is not compatible with the Langevin model, 
which assumes the moments are free to rotate. In this 
respect, weakly uncompensated antiferromagnetic parti- 
cles strongly differ from ferro- or ferrimagnetic particles. 
In the latter, the anisotropy energy is easily overcome by 
the Zeeman energy and the Langevin model is a good 
approximation, although deviations from it have been 
observed at low temperature and interpreted using a full 
Boltzmann calculation of the magnetization (see for in- 
stance Ref. f26l l ) . In randomly oriented antiferromagnetic 
particles, a calculation of the total magnetization along 
the same lines is much more difficult because it involves 
four variables, i.e. two positional angles for the each of 
the two sublattices. Therefore, as a low field approxima- 
tion, we interpret the contribution of the uncompensated 
moments in the superparamagnetic regime using the dy- 
namic "random magnetic orientations" model, assuming 
that they fluctuate between two antiparallel directions 
along the antiferromagnetic axes. The uncompensated 
moment, for a given particle volume and a given orienta- 
tion, is then described by a hyperbolic tangent function 
similar to the expression for the hyperfine field (Eqn. 12), 
and the measured moment is obtained by integration over 
the random directions of the antiferromagnetic axis: 



3 (T, V)G 



./x nc (T, V)H. 



k B T 



with the function G(x) defined by: 



G(x) 



1 



d6 sin 6 cos 6 tanh(a; cos ( 



(16) 



(17) 



The G(x) function is plotted in the inset of Fig.fLO] to- 
gether with the Langevin function L(x); G(x) saturates 
at the value 1/2, whereas the Langevin law saturates at 
unity. The isothermal magnetization curves have been 
fitted with this model of dynamic "random magnetic ori- 
entation" , according to the law: 



M{T,H) 
1 



< V > 



Xaf(T)H 

V ma x 



ri nc (T,V)f(V)G 



u. nc (T,V)H 
k B T 



]dV, (18) 



where f(V) is the lognormal volume distribution func- 
tion. We chose to describe the uncompensated moments 
by their mean value: fJ, nc (T, V) = mou(T)y/N, where 
u(T) accounts for the thermal variation of a Fe 3+ mo- 
ment. In Fig.|H^ are represented the low field parts of the 
isothermal magnetization curves in the particles with 982 



Fe atoms per core, at 35 K and 80 K, together with the fits 
of the data using the G(x) function (solid line) and the 
Langevin function (dashed line) with the same volume 
distribution. It is clear that the fit using the "random 
magnetic orientation" model gives a better account of 
the M(H) curvature at low fields than does the Langevin 
fit, the improvement being more pronounced at low tem- 
perature. The complete field and temperature variation 
of the magnetization in the particles with 982 Fe atoms 
per core is shown in Fig. 
the fit using the 
(Eqnjl|). 
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the solid lines representing 
"random magnetic orientation" model 
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FIG. 11. Magnetization curves in the artificial ferritin sam- 
ple with a mean Fe loading of 982 atoms per core. The solid 
lines are obtained by using the "random magnetic orienta- 
tion" model (see text, Egn.^Ufy) The inset shows the thermal 
variation of xaf in the same sample. 



These fits are very satisfactory in the whole field range, 
and the values for u are about 0.6 (close to those obtained 
in the Mossbauer spectra at 90 K at the lowest fields) and 
slowly decrease as temperature increases. We emphasize 
that the description of the thermal and field variation 
of the uncompensated magnetization in the superparam- 
agnetic regime using the G(x) function lies on stronger 
physical grounds than that using a Langevin function. 
The upper bound of the field range where this approach is 
valid, as discussed in section VI, is approximately Ha/cx, 
i.e. it is close to the maximum field (5-6 T) used so far in 
magnetization studies of ferritin. Above this field, the un- 



compensated magnetization starts to decrease and van- 
ishes when the field reaches the spin-flop threshold, as 
the uncompensated moments fluctuate along a direction 
close to perpendicular to the applied field. The G(x) 
function has been evoked by Neel jl6],^5| in his studies 
of antiferromagnetic particles: he applied it to calculate 
the thermoremanent magnetization at the blocking tem- 
perature and stated that it was the correct description of 
the uncompensated magnetization in the limit of strong 
anisotropy. This fact has also been noted by Bean |29] 
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for the magnetization of ferromagnetic fine powders. We 
showed here that, in weakly uncompensated antiferro- 
magnetic particles (i.e. with a degree of uncompensation 
of a few percent), the G{x) function yields a correct de- 
scription of the magnetization even for finite and rather 
small values of the anisotropy energy, like those found 
in ferric oxides, for fields up to 5-6 T. This is finally a 
consequence of the very large value of the exchange field, 
which yields a large spin-flop field. 

The thermal variation of the mean uncompensated mo- 
ment derived from our fits of the M(H) curves presents 
a small but significant decrease as temperature increases 
from 35 K to 250 K, for all particle sizes. Using an anti- 
ferromagnetic magnon law: /i„ c (T) = /i nc (0)(l — otT 2 ), 
the experimental /i nc (T) curves can be extrapolated to 
zero in order to yield an estimation of the Neel temper- 
ature. We find that Tm is close to 500 K, and that it is 
practically the same for all particle sets. This absence 
of dependence of TV on the mean particle size is prob- 
ably due to the fact that the distance between "active" 
reticular planes, which determines the strength of the 
superexchange interaction, is independent of the parti- 
cle size. The Tn value we obtain is roughly two times 
larger than early estimations |]||], but is ingood agree- 
ment with more recent determinations |^,|9|,|l0[ in natural 
ferritin. 

The antiferromagnetic susceptibility \af{T), shown in 
the inset of Fig.[ll] for the particles with 982 Fe atoms 
per core, decreases with increasing temperature, whereas 
the bulk powder antiferromagnetic susceptibility is ex- 
pected to increase on heating. A mechanism for the 
thermal decrease of xaf has been proposed by Neel [^5) , 
by considering the reduced exchange field of the super- 
ficial Fe layers, but it is difficult to say whether it is 
responsible for the thermal decrease of xaf we observe 
in the ferritin particles. The T=0 value of xaf, de- 
rived from the magnetization curves at 2.5 K (not shown 
here), amounts to a few 10 emu/cm 3 for all particle 
sets. In the simple molecular field model, the expression 
for the powder T=0 antiferromagnetic susceptibility is: 
Xaf(T = 0) = |x_l = fir' Using the estimated ex- 
change field value -£/e=320T, we find that, for each par- 
ticle set, the calculated xaf(T = 0) is smaller than the 
experimental value by a factor of about 3. Neel showed 
that the susceptibility of antiferromagnetic particles hav- 
ing an even number of "active" reticular planes (i.e. no 
net magnetic moment) can be enhanced by a factor 2 
or more with respect to the bulk susceptibility, and that 
this is a surface effect visible in very small particles |30) . 
He called this behaviour "superantiferromagnetism" . We 
think that this mechanism is effective for such particles 
even in the presence of a small uncompensated moment, 
and that it is responsible for the observed large value 
of the antiferromagnetic susceptibility. The superantifer- 
romagnetic behaviour is confirmed by high field (30 T) 
magnetisation measurements we performed at 2.5 K in 
natural ferritin, which will be the subject of a future 
publication pl[ . 



At this stage, we would like to comment on the ob- 
servability of effects linked with surface atoms in antifer- 
romagnetic particles. Due to the small size of our parti- 
cles (4-6 nm), an important fraction of Fe atoms (around 
40%) lie at the surface, and their properties (anisotropy, 
exchange field) are likely to be somehow different from 
those of the core atoms. For instance, in ferro- or fer- 
rimagnetic systems, the observed lack of full alignment 
of the magnetic moments with the applied field is char- 
acteristic of nanoparticles and is generally attributed to 
the different behaviour of surface atoms jl8|,[32| ■ In anti- 
ferromagnetic particles in moderate magnetic fields, the 
situation is different because the dominant configuration 
is an almost random orientation of the moments with 
respect to the field. Then, effects linked with surface 
moments are much more difficult to evidence experimen- 
tally and to distinguish from those due to core moments 
(see however Ref. [p3|). So, although the properties of 
surface atoms may be different from those of the core 
atoms in our artificial ferritin particles, we could not de- 
tect their effects in our experimental data, except for the 
superantiferromagnetic behaviour mentioned above. 



VIII. CONCLUSION 

We performed 57 Fe Mossbauer absorption spec- 
troscopy and magnetization measurements in antiferro- 
magnetic artificial ferritin particles, with mean Fe load- 
ings ranging from 400 to 2500 atoms, in the tempera- 
ture range 2.5K-250K and with magnetic fields up to 
7T. In zero or very low field, the dynamics of the sub- 
lattice magnetization of the ferritin particles obeys clas- 
sical superparamagnetic relaxation in the thermal acti- 
vation regime. The value of the anisotropy energy per 
unit volume K could be determined from the FC-ZFC 
susceptibility curves and from the thermal variation of 
the Mossbauer superparamagnetic fraction. The values 
obtained by both techniques are in good agreement and 
are in the range 3 — 6 x 10 5 ergs/cm 3 , typical for ferric 
oxides or hydroxides, with a tendency for K to increase 
as the particle size decreases. By comparing the data ob- 
tained by Mossbauer spectroscopy and by magnetometry 
in moderate magnetic fields up to 7T, we propose a new 
interpretation of the field and temperature behaviour of 
the uncompensated magnetization in the superparamag- 
netic regime. We show that it is better described by a 
"random magnetic orientation" model, where the uncom- 
pensated moments fluctuate along the antiferromagnetic 
axis, than by the usually assumed Langevin law. This 
behaviour is described by the universal function: 

1 r 

G(x) = - d9 sin 9 cos 6»tanh(ic cos 9). (19) 

2 Jo 

To our knowledge, this law has never been applied to 
the interpretation of data in antiferromagnetic particles, 
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although Neel had noticed that it should be the cor- 
rect description of the thermal and field variation of 
the uncompensated magnetization in the limit of large 
anisotropy |25f| . We think this non-Langevin behaviour 
actually holds for usual values of the anisotropy energy, 
and it should be valid not only in ferritin, but also in 
other weakly uncompensated antiferromagnetic particles, 
at least in moderate magnetic fields. As the field is in- 
creased towards the spin-flop field (10-20 T) and above, 
the "random magnetic orientation" picture breaks down, 
as the magnetic moments progressively reorient perpen- 
dicular to the applied field and the superparamagnetic 
uncompensated magnetization vanishes. 

We could also estimate the Neel temperature in our 
artificial ferritin samples by following the decrease of the 
mean uncompensated moment in the temperature range 
35-250 K; we find that Tn is essentially independent of 
the mean particle size, and that it is close to 500 K, in 
agreement with other recent estimations. 

The authors are grateful to Dr E. Vincent and Dr G. 
LeBras, from the Service de Physique de l'Etat Condense 
(CEA Saclay), for their help with the SQUID measure- 
ments. 
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